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Abstract This is a brief review of few relevant topics on tunneling of composite particles and how 
the coupling to intrinsic and external degrees of freedom affects tunneling probabilities. I discuss the 
phenomena of resonant tunneling, different barriers seen by subsystems, damping of resonant tunneling 
by level bunching and continuum effects due to particle dissociation. 
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1 Particles moving in mysterious ways 

History - The tunnel effect is one of the most subtle phenomenon explained by quantum mechanics, 
responsible e.g., for the existence of stars and ultimately for the existence of life [I]. The first application 
of particle incursion into classically forbidden regions was done in nuclear physics. Nuclei, such as 210 Po, 
emit o-particles by tunneling through the Coulomb barrier. This process lacked physical explanation 
until George Gamow used the tunneling theory to calculate o-emission half-lives El Eg. Tunneling 
is now a well known physical process that has been incorporated in our everyday lives due to the 
increasing miniaturization in electronics such as microchips. There have been important Nobel Prizes 
related to the straightforward use of this effect in industry, such as the tunnel diode [4] or the scanning 
tunneling microscope [5]. 


Resonant tunneling - Resonant tunneling is a particular kind of tunneling effect, frequently applied to 
miniaturization such as the resonant diode tunneling device. In its simplest form, resonant tunneling 
occurs when a quantum level in one side of a barrier has an energy match with a level on the other side of 
it. If this occurs in a dynamical situation, tunneling is enhanced. In certain situations, the transmission 
probability is equal to one and the barrier is completely transparent for particle transmission. In the 
resonant diode tunneling device, two semiconductor layers sandwich another creating a double-humped 
barrier which enables the existence of quantum levels within. On both sides of the barrier, electrons fill 
a conducting band on two outside semiconductors. A potential difference between these outer layers 
allows one conducting band to rise in energy. Resonant tunneling leaking occurs to the levels inside 
the double-hump barrier [4]. The net effect is the appearance of a current which can be fine tuned. 
Resonant tunneling devices are compact and allow a fast response because the tunneling between the 
thin double-humped barrier is a very fast process. 
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Fig. 1 Left: Tunneling of a composite particle through a barrier is often equivalent to the tunneling of a single 
particle through a double-humped barrier. Right: A Feshbach resonance appears when a preferred tunneling is 
induced by the coupling of a reaction channel to one or more intrinsic degrees of freedom. 



Fig. 2 Left: Tunneling of a composite particle through a single barrier. Right, (a) and (b): The effective 
potential for the center of mass motion of the system. 


Composite particles - Resonant tunneling is not constrained to a particle tunneling though a double- 
humped barrier. An equivalent process occurs when a composite particle tunnels through a single 
barrier, if each subsystem composing the particle is allowed to tunnel independently or when the 
interaction reveals a priori unknown energy states. Then the process is equivalent to the tunneling 
of one single particle though a double-humped barrier, and resonant tunneling proceeds through the 
appearance of pseudo, quasi-bound, states (see Figure IT] left). This effect occurs frequently in atomic, 
molecular and nuclear systems. A good example is the fusion of loosely-bound nuclei [6]. It also relates 
to Feshbach resonances when a system uncovers a state in a (closed) channel which is not the same 
(open) channel where it sits in. This happens because the coupling with at least one internal degree 
of freedom helps the reaction to proceed via resonant tunneling (see Figure [l] right). 

In the next Section I give several examples of how intrinsic degrees of freedom manifest in subtle 
ways, leading to tunneling enhancement or suppression. As the subject is very vast and general, I will 
concentrate on examples that I had the opportunity to work with during the last decades. They range 
from particle to molecular physics. Some common characteristics are easy to understand, others not so 
much. An example of the last case if given for diffusion and dissociation of molecules in optical lattices. 


2 Tunneling of composite particles 

Most theoretical problems involving tunneling are not amenable to analytical solutions. There are 
several computational methods available for them. Below I will just mention a very well-known method, 
used to generate some of the calculations displayed in the figures of this review. 


Stalking microscopic particles - It is interesting to study the time-dependence of tunneling, although 
it is prohibitive for most cases of interest due to the long tunneling times, e.g., in nuclear <a-decay 
processes. The time evolution of a wave function on a space lattice can obtained by solving the 
Schrodinger equation by a finite difference method. The wave function \P(t + At) can be calculated 
from the wave function at time £, ^(t), by applying the unitary time evolution operator, U, i.e., 
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Fig. 3 Left: Transmission probabilities for a two-body system through a rectangular barrier m as a function 
of the ratio of the center of mass energy and the barrier height Vo- The solid line is the transmission probability 
in which the system remains in the ground state and the dashed line is the for the case the system transits 
to an excited state. Right: Transmission probability of a molecule through a thin barrier as a function of the 
barrier width a times the molecule momentum. The dashed curve is the transmission probability when the 
molecule is in an excited state after the transmission na- 




Fig. 4 Left: A schematic rectangular barrier simulates the proton-nucleus interaction for a deuteron-nucleus 
fusion reaction. The solid curve is the effective potential for the center of mass motion. Right: Transmission 
probability for the deuteron, with (solid) and without (dashed) consideration of coupling to intrinsic modes 


T(t + At) = U(At 1 t)T(t) = exp [—iHAt/h]\P(t), where H is the system Hamiltonian. For a small time 
step At , one can use the unitary operator approximation, valid to order ( At ) 2 , 


U(t + At) 


1 + ( At/2ih)H(t) 

1 - ( At/2ih)H(t) ’ 


(i) 


Other approximations for exp [—iHAt/h] such as the Numerov algorithm, etc., can also be used. Any of 
these numerical procedures requires carrying out matrix multiplications and inversions at each iteration 
and the number of operations grows very fast with the number of points in the coordinate lattice 00. 
It also increases with the number of particles in the system at hand, becoming numerically demanding 
for large fermionic or bosonic systems [9]. But the most difficult problem is related to the existence of 
times scales which differ by a large amount. For example, during alpha decay the decay time is often 
much larger (e.g., 10 30 ) than the intrinsic time dependence of the alpha particle state in the nucleus. 
This makes it impossible to treat the problem in a time-dependent fashion. But there are ongoing 
theoretical efforts to handle the time-dependent problem with widely different time scales m nn US¬ 
As an application of the method described above, I will present below results for the tunneling of 
a two-body system through a simple rectangular barrier. 


Transmission of composite particles - Assume a system of two bound particles with an attractive 
infinite square-well interaction of range d, being transmitted through a rectangular barrier of width 
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a and height Vo. By choosing the center of mass x = (x\ + xf)/2 and relative y = x\ — + d/2 

coordinates, it is straightforward to show that the center of mass wave equation has an effective 
double-humped potential W(x, y ). Figure [ 2 ] shows that for various values of the relative distance y , the 
effective potential W is in fact a double-humped barrier 0H2]. This barrier can hold energy states 
that yields resonant tunneling if they are equal to the center of mass energy, as shown in Figure [3] for 
the dependence of the ratio of the bombarding energy and the barrier height SHU]. 

The simple example mentioned above can be readily applied to a very practical case; that of a 
bound molecule transmitting through a barrier. This is illustrated below. 

Transmission of molecules - Another way to enhance tunneling is allowing the system to make a 
transition from an excited state to a lower energy state. This has been shown analytically in Ref. m 
for a two-spin system tunneling though a barrier. The transmission probability as a function of the 
barrier width also shows signs of resonant behavior, as observed in calculations for transmission of 
molecules though barriers m- Figure [ 3 ] shows the transmission probabilities for a two-body system 
through a rectangular barrier [121 as a function of the ratio of the center of mass energy and the barrier 
height Vq- The solid line is the transmission probability when the system remains in the ground state 
and the dashed line is when the system transits to an excited state. The right hand side of Figure [3] 
shows the transmission probability of a molecule through a thin barrier as a function of the barrier 
width a times the molecule momentum. The dashed curve is the transmission probability when the 
molecule is in an excited state after the transmission m- By moving to its ground state the excited 
molecule increases its chance of transmission. On the opposite trend, the opening of the higher channels 
decreases the probability of transmission. 

The two-body tunneling problem also encounters abundant applications in nuclear physics, involv¬ 
ing cluster-like objects. The simplest one is the deuteron, composed of a proton and a neutron forming 
a relatively loose-bound system (Eb = —2.224 MeV). Other common systems include those containing 
alpha particles. Tunneling involves sudden accelerations leading to gamma-emission in alpha tunneling, 
as described below. 

Fusion of loosely bound nuclei - Loosely-bound systems also display interesting tunneling properties 
because they can end up in separate parts by reflection or transmission through a barrier. If the 
separate parts see a different barrier then the effective barrier for the center of mass tunneling can 
also be very different than the original one. A good example is the fusion of the deuteron with another 
nucleus in which case the proton-nucleus potential has a Coulomb repulsion part whereas the neutron- 
nucleus potential does not. As shown in Ref. [14] the net effect is an enhancement of the transmission 
probability through the barrier (see Figure®. Fusion of exotic, neutron rich, nuclei due to particle 
transfer, or tunneling, is now a very active area of research [16.; 17; 18; 19; [20; [2TJ [22; [23] (for a review, 
see Ref. [24] 1. 

An intriguing question is how much time a particle takes to tunnel a barrier [25]. A Bremsstrahlung 
measurement following a -decay m apparently paved the way to determine tunneling times by inspect¬ 
ing the interference pattern in the radiation spectra. A time-dependent description of the emission pro¬ 
cess with a pre-formed <a-particle has shown that a model with a free, uncorrelated, <a-particle within 
the nucleus is incompatible with the observations for this process E3- The spectrum of Ref. [26] was 
also explained in a more traditional theoretical method ESI Ea EDI. However, the question related to 
tunneling times is still fascinating, in particular when it comes to the possibility of time travel and the 
existence of time warps and worm holes m- 

The emission of radiation during tunneling is also manifest in virtual processes, a fascinating dis¬ 
cussion of which is presented below. 

The quantum Miinchhausen effect - Additional interactions during tunneling can lead to other subtle 
effects, e.g., as proposed in Ref. [321 • Since barrier transmission is proportional to an exponentially 
decreasing dependence with the particle energy, any small effect can change it dramatically. The so- 
called quantum MunchhauserQ effect assumes that a virtual photon emitted in one side of the barrier 
and caught on the other side can enhance barrier penetrabilities (Figure |5|. This is allowable in 

1 Baron Miinchhausen is supposedly the greatest liar ever depicted in fiction literature and in movies. He 
once explained how he could save himself from sinking in quicksand by pulling his own hair up [33] . 
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due to couplings to real and virtual photons. 




He + H + H + b (a.U.) 


Fig. 6 Left: Energy levels of an electron as a function of the distance between a proton and a hydrogen atom. 
Right: The exchange probability as a function of the impact parameter for two proton kinetic energies. 


quantum mechanics because a particle has an extended wave nature. The particle kinetic energy can 
momentarily increase (within bounds set by Heisenberg’s uncertainty principle) and tunneling can be 
enhanced. A non-relativistic reduction of the particle-radiation interaction yields [32 

H#(r) + J n(r,T f -,E)&(T f )dr f = E&(r), where i7(r, r'; E) = M(r, r'; E) + ir(r, r'; E). (2) 

In this expression, M includes the self-energy, mass renormalization and virtual photons, while r 
includes the effects of radiative decay width and real photons. This effect has been studied further in 
Ref. [34] and it was shown to be very small indeed for tunneling of stable nuclei. But an account of 
this effect for loosely-bound nuclei (and molecules) has not been explored at depth. 

In atomic systems, electron tunneling is sometimes responsible for chemical bonding. But subtle 
effects are also present in dynamical situations, in which electron tunneling can be the cause of increased 
stopping of protons on a hydrogen target at extremely low energies, as discussed next. 

Electron hoping induces stopping - As discussed above, tunneling can be strongly affected by matching 
energy levels, the appearance of pseudo-levels during the tunneling process (Feshbach resonances), and 
by energy transfer to the relative motion by coupling to the internal degrees of freedom. Other effects 
appear when many energy levels are involved. Perhaps the best and clearer examples appear in atomic 
physics. Stopping power, or the energy loss of a projectile per unit length, S = — dE/dx , is an important 
quantity for experiments at low bombarding energies. The energy at which a nuclear reaction occurs 
depends on the average energy loss of the projectile in an atomic target. Correcting for the stopping 
power is thus imperative for low energy experiments. At extremely low projectile energies, only a few 
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Li+ He + H + 




Fig. 7 Left: Adiabatic energies (1 a.u. of energy = 27.2 eV, 1 a.u. of length = 0.53 A) for the electronic orbitals 
of the (H-He) + system as a function of the internuclear separation. As the atoms approach each other slowly, 
curves of same symmetry repel each other. A transition between states s and s’ can occur in a slow collision. 
In a fast collision a diabatic transition, with the states crossing each other, will occur. This is shown in the 
inset. Right: Probability of charge exchange in the collision p + 4 He showing the damped resonant behavior as 
a function of the impact parameter and for proton energy E v — 10 keV. 


phenomena are energetically viable such as the excitation of vibrations (plasmons) in the medium for 
solid materials, Rutherford scattering and straggling, or charge exchange when one or more electrons 
are transferred to the projectile. 

Stoping by charge exchange in p + H and p + D reactions has been studied in Ref. [35] , Because 
of the low bombarding energy, the electrons adjust themselves quickly and one can use the adiabatic 
approximation. The potential felt by the electron resembles a symmetric Coulomb barrier around 
the two atomic centers. As a function of the relative nuclear distance, the electronic energy levels 
smoothly undergo from hydrogen to helium levels (left side of Figure [6]), passing by molecular orbitals 
at intermediate distances. Expanding the electronic wavefunction as a time admixture of such orbitals, 
a set of coupled-channel equations is obtained with help of the Hellmann-Feynman relation. As the 
bombarding energy decreases only the two lower levels, Isa and 2 pa play a relevant role. The electron 
jumps back and forth between the two atoms by tunneling through the symmetric Coulomb barrier. 
The right hand side of Figure [6] shows the electron exchange (tunneling) probability as a function of the 
collision impact parameter in atomic units and for two bombarding energies: (a) 10 KeV, and (b) 100 
eV. One notices that at higher energies the tunneling probability increases as the distance decreases 
mainly due to the decrease of the barrier height and the resonant tunneling effect. When the energy 
is very low, as seen in the lower panel, the electron jumps frenetically back and forth between the two 
atoms due to resonant tunneling. At very low energies the exchange probability as a function of impact 
parameter is given in terms of the simple formula [35] 

Pexch = 7} + 7} cos 11 J [E 2p (t) - E u (t )] dt | . (3) 

The dashed curve in Figure [6] proves that the prediction based on this equation holds very well, 
specially at very low energies. It is expected to be exact at extremely low energies. The probability 
minima occur for impact parameters satisfying the relation [E 2 P (t) — Ei s (t)] dt = 2 tt h(n+ 1/2), a 
relation more than familiar to all of us. It is the interference between the Isa and the 2pa states that 
induces oscillations in the exchange probability. 

The case discussed above involves tunneling of only one electron. A more interesting and rather 
complex situation arises when two electrons are present, as in the case of stopping of protons on a 
helium gas, as discussed below. 
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Damping of resonant tunneling - The case of p + He is more interesting, because of the two active 
electrons. One can expand the electron wave functions in terms of the two-center basis and obtain a 
set of time-dependent Hartree-Fock equations [36] . In contrast to the p + H or p + D case the much 
larger number of levels make a difference. Due to level repulsion, the levels tend to bunch as seen in 
the left panel of Figure [7] The transition between the lowest levels shows an oscillating behavior due 
to resonant tunneling between one electron in the ground state of hydrogen and the first excited state 
of helium. But the oscillations are damped. This is due to the interference between the low-lying states 
and a bunch of states of average energy (E) and width 2T, as shown on the left side of Figure [7] The 
damping can understood with the Landau-Zener theory for level crossing. At the crossing there is a 
probability of an adiabatic transition where (1 — P eX ch ) given by the Landau-Zener formula 


Pexch = exp 


2ttH 2 , 


_kvd(E s — E s >)/dR_ 


( 4 ) 


where v is the collision velocity, H ss > is the off-diagonal matrix element connecting states s and s' and 
the term in the denominator is the derivative of the energy difference with respect to the separation 
distance between the two nuclei. The interference with the neighboring states introduces a damping in 
the charge exchange probability, 


Pexch — COS 



2irr 2 b\ 

hv (E)) 5 


( 5 ) 


where ( E ) « 1 a.u. is the average separation energy between the OE level and the bunch of higher- 
energy levels shown in Figure [7] The exponential damping factor agrees with the numerical calculations 
if one uses r = 5 eV. The exchange probability for the electron drops sharply to zero when the energy 
of the proton is smaller than 18.7 eV. This is easy to understand because this is the energy necessary 
for the transition Is 2 (-'-So) — 1s2s( 3 S) in He. Only when this transition occurs, resonant tunneling 
will be possible between the ground state in hydrogen and the excited state in helium. 

Finally, we consider what happens in tunneling of composite objects through multiple barriers. As 
an illustrative example, we describe below the case of diffusion and dissociation of molecules in optical 
lattices. Despite the theoretical simplicity of the system under consideration, the complexity is much 
larger than in the previous examples and some open questions remain. 


Diffusion in optical lattices - Using standing wave patterns of reflected laser beams optical lattices can 
be built to study the diffusion of ultracold atoms and molecules [371131ES ED]. Fundamental aspects 
of diffusion of composite objects, such as molecules, subject to transforming transitions from bound 
to continuum states can be studied in these optical lattices [41]. In Ref. [41] the confining potential 
described by the optical lattice potential has been assumed to be a sine function with a periodicity 
D, equal to half the wavelength of the interfering laser. The Hamiltonian for the molecule is given by 
H = Ti + T 2 + Vi + V 2 + v, where Ti is the kinetic energy of atom i, Vi is the periodic potential of the 
lattice and v is the interaction potential between the atoms for a diatomic molecule. A loosely bound, 
Rydberg-like, system of molecules diffusing through the lattice molecule was studied in Ref. m- 

The molecules tend to dissociate as they tunnel through the barrier when their binding energy 
decreases. At very small binding energies, the trend is reversed and the molecule tends to remain 
intact. As shown in Figure [8] loosely-bound molecules are resilient to breakup: even at the smallest 
binding value, the breakup probability is not one. This result is not unusual, as other tunneling systems 
exhibit similar behavior. This has been proven for the familiar example of tunneling of Cooper pairs 
[42l . The pair does not usually dissociate as it tunnels through a barrier. However composite particles 
may dissociate as the particles tunnel through multiple barriers, as in the case of a lattice potential. 
It takes time for the molecules to dissociate and additional time for the atoms to diffuse after the 
dissociation. Therefore, the molecules initially confined within one pocket of the lattice will tunnel and 
diffuse away from the initial position at a speed higher than that of the atoms that are created in the 
dissociation of the molecules. 

In Ref. [41] a diffusion equation for molecules in an optical lattice in terms of the diffusion velocity 
of the wave packet was given by a(t) = yj (r 2 (t)), where (r 2 (t)) is the expectation value of the square 
of the position of the molecule calculated from the solution of the Schrodinger equation. The numerical 











Fig. 8 Left: Lattice diffusion of molecules. The hatched histograms are the relative probabilities of finding 
a molecule in its ground state at a given position along the lattice. The solid histograms give the relative 
probabilities of finding individual atoms after the dissociation. The molecules are initially confined in the site 
at the origin. Right: Time dependence of the spreading width of bound molecules, <jm(£), shown by the solid 
line, and of dissociated atoms, cta(£), shown by the dotted line. The dashed curve is a fit of the asymptotic 
time dependence with the analytical formula, Eq. <§• 


results show that after a transient time, and for a strongly bound molecule of mass ra, the approximate 
result holds, 


a(t) = C 


h 2 t 

mbcjQ 


1/2 


( 6 ) 


where C and b are constants depending on the parameters defining the lattice and ao = a(t = 0). b 
has the meaning of a “viscosity” generated by the wiggling potential dependence of the lattice. The 
linear behavior in Eq. ^ is shown by the dashed line in Figure [8] (right). Also shown are the numerical 
results for the same quantity (solid line), and the diffusion velocity of individual atoms in the lattice 
(dotted line). 

A diffusion coefficient can be defined by D = (da 2 / dt)/2. Using Eq. ^ one gets 




( 7 ) 


This result deviates from the classical Einstein diffusion constant, D = kT/b. The effective “tem¬ 
perature” for this system is the kinetic energy of the initial state, T ~ Ti 2 /ma 2 (t = 0). The dependence 
of the temperature, or diffusion coefficient, on the initial kinetic energy is a result which needs to be 
verified experimentally. The molecules are assumed not to interact. Therefore these results might not 
hold when molecules interact in a condensate. The time dependence of the spreading width of the 
dissociated atoms within the lattice is not well fitted by any power-law dependence on time, as shown 
in Figure [8] The reason for this behavior is not yet understood. 


3 Conclusions 

In this short review I have presented an overview of tunneling of composite objects. I started by giving 
an example of how one would tackle the problem computationally, followed by an application of the 
method to the tunneling through a rectangular barrier. It is evident that coupling to open and closed 
channels lead to a clear resonant behavior visible in the tunneling probabilities. 

I went on to discuss a few applications of tunneling of composite objects to real situations, such as 
molecules passing through thin membranes, cluster-like nuclei tunneling during a nuclear fusion process, 
and the emission of radiation during tunneling as, e.g., observed during alpha-decay. Other examples 
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were given in atomic physics where tunneling can enhance or suppress the stopping of particles at 
very low energies. This important process leads to corrections for the energy loss of ions in reactions of 
astrophysical interest. Finally, the more difficult theoretical treatment of tunneling of composite objects 
through multiple barriers has been illustrated with the example of molecules diffusing in optical lattices. 

Although a (large) number of theoretical works have studied tunneling phenomena in various situ¬ 
ations, quantum tunneling of a composite particle in which the particle itself has an internal structure, 
has yet to be fully clarified. 
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